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Abstract 

For the Landau problem with a rotating magnetic field and a potential in the (changing) direction 
of the field, we derive a general factorization of the time evolution operator that includes the adiabatic 
factorization as a special case. The potential is assumed to be of a general form and it can correspond to 
nonlinear Heisenberg equations of motion. The rotation operator associated with the solid angle Berry 
phase is used to transform the problem to a rotating reference frame that follows the direction of the 
magnetic field. In the rotating reference frame, we derive a natural factorization of the time evolution 
operator by recognizing the crucial role played by a gauge transformation. The major complexity of 
the problem arises from the coupling between motion in the direction of the magnetic field and motion 
perpendicular to the field. In the factorization, this complexity is consolidated into a single operator 
that approaches the identity operator when the potential confines the particle sufficiently close to a plane 
perpendicular to the magnetic field. The structure of this operator is clarified by deriving an expression 
for its generating Hamiltonian. The adiabatic limit and non-adiabatic effects follow as consequences of 
the general factorization which are clarified with the magnetic translation concept. 



1 Introduction 



The study of time-dependent quantum systems has intimate connections with the geometric phase concept 
[TJ[3J[3] which has many applications in physics. When the Hamiltonian is time-dependent, the time evolution 
operator is often nontrivial, i.e., U{t) ^ exp(— i H{T)dT), and clariiying its structure is of significance in 
understanding the dynamics of the system. 

One may study U (t) corresponding to a time-dependent Hamiltonian by factorizing it into several op- 
erators, each of them is simpler at least in some respects than U (t) itself. A well-known example of this 
is found in the proof of the quantum adiabatic theorem, as presented in standard texts such as Messiah 
[4] . There the time evolution of a system in a changing environment is constructed as the product of three 
operators, U = GDUe, where G is a path dependent geometric operator that brings an initial eigenstate to 
an instantaneous eigenstate of the Hamiltonian, D is a dynamical operator that only contribute dynamical 
phase factors to these eigenstates and Ue approaches the identity operator in the adiabatic limit e — >■ 0. The 
parameter e determines how fast the Hamiltonian _ff [R(s)] — H[EL{et)] changes with time t, for a given map 
from (s €)[0, 1] to a path in the parameter space M to which R belongs. For any such fixed map from [0, 1] 
to M, 1/e provides a time scale which is the total time it takes for R(ei) to travel through the given path in 
M. Another relevant time scale is provided by h/Eg, where Eg is the minimum energy gap. The adiabatic 
limit corresponds to {h/Eg)/{l/e) — ^ 0, which is made use of in the proof of the adiabatic theorem [3J [S]. 
The path-dependent geometric operator is (assuming the Hamiltonian has non-degenerate eigenstates for all 
times) 

G(R(ei)) = IV'm(R(rf))>(V'™(R(0))|, (1) 

where |'0TO(R(e<))) is the instantaneous eigenstate oi H\R,{et)] that satisfies 

(V'™(R(ei))l^™(R(eO)> =0. (2) 

The expressions for the operators D and also involve the use of eigenstates of the Hamiltonian. 

This method of factorizing U using instantaneous eigenstates has certain limitations when dealing with 
degenerate (including infinitely degenerate) energy eigenstates. For instance, in various Landau systems 
involving a charged particle in time-dependent electromagnetic fields, the instantaneous energy levels can 
be highly degenerate. Then there is no known general method for obtaining useful information on G (which 
contains information on non-Abelian Berry phase) or using instantaneous eigenstates. However, in a 
specific problem where the Hamiltonian is given, one may use the algebraic structure of the Hamiltonian 
without referring to individual eigenstates to directly construct a factorization of U that can then be applied 
to any representation and the associated eigenstates. From this perspective, there seems to be more problems 
that can be explored. 

This change in perspective also allows us to seek useful factorizations of U not limited by the specific form 
U — GDUe, as long as the factorization can help clarify the structure of the total time evolution operator U. 



2 The problem and general considerations 

The Landau problem is of significance in many areas in physics and its variations (see, for instance, [BJ [71 ISl IS] ) 
have often been discussed. In this paper our purpose is to study a charged particle in a rotating magnetic 
field and a confining potential in the direction of the magnetic field. The Hamiltonian is 

H=^{p- eA(r, t)f + V{r ■ n{et) - L), (3) 

where 

A(r,t) = -Bn(et) X r. (4) 
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The confimng potential V(r ■ n{et) — L) is in the (changing) direction of the magnetic field Bn{et) and has 
equilibrium position at the plane r • n(et) — i = which is perpendicular to n{et). The distance between this 
plane and the origin of the coordinate system is L. This can be seen as an extension of the usual Landau 
problem where n is in a fixed direction. 

Here, y(r-n(et) — i) is assumed to be of a general form. It can be a harmonic oscillator potential or other 
types of potentials that in general correspond to nonlinear Heisenberg equations of motion. The complexity 
of the problem is mainly caused by the coupling between motion in the magnetic field direction and motion 
perpendicular to it. The question we ask is whether the time evolution operator has a factorization that 
reduces to a simple form when the confinement is strong enough so that the particle can be seen as staying 
very close to the rotating plane. There are three time scales involved in the problem: 

Ti = l/e, T2=27r/io, Tg = ^/A, (5) 

where uj — \e\B/m is the cyclotron frequency and A is the minimum energy gap determined by the Hamil- 
tonian 

±-(p.nf + V{v-niet)-L). (6) 

The general factorization we derive later does have a simple form when the confinement is strong in the sense 
that Ti and <^ This factorization is from a different perspective than a factorization that results 

in a simple behavior when e is small (i.e., T3 <^ Ti and T2 Ti). Here we point out that this different 
perspective seems to be more natural because, in this problem, it seems more appropriate to consider the 
adiabatic limit to be a special case of the strong confinement limit. In other words, the adiabatic limit here 
should mean 

Tz/Ti-^O, T3/T2^0. (7) 

The conditions T2/T1 — )■ and T3/T1 — are weaker than the above and do not lead to a simple behavior 
of the time evolution operator, since if T2 is comparable to T3, the coupling effect between motion in the 
direction of n and motion perpendicular to it can accumulate over the time interval [0, Ti] even if e is very 
small. This was pointed out in a special case [7] where we studied the factorization of U in the adiabatic 
limit only and where we took V to be the harmonic oscillator potential V = ■^k(r ■ n{et) — L) . The 
special case, however, misses out on the more general scenario of strong confinement where U already has a 
simple factorization. Furthermore, the adiabatic factorization contains no information on effects due to finite 
confinement strength or that due to non- vanishing rotation speed. Using an improved method, in which we 
recognize the key role played by a gauge transformation, our purpose here is to derive a factorization of U 
that applies for all confining potentials and an arbitrary time variation of n(ei). 

The electric field E(i) = — ^ corresponding to the vector potential ^Bn{et) x r is a notable point 
although the potential is quite commonly used in the literature. (Here we assume that the magnetic field 
induced by the change of E(i) and so on are negligible.) The subtlety lies in that the vector potential is 
time dependent, unlike a vector potential that corresponds to a static magnetic field. Imagine we choose a 
different vector potential through a time dependent gauge transformation, then we have the same rotating 
magnetic field but a different electric field. Besides considerations based on symmetry and simplicity, there is 
no physical reason why we prefer this vector potential to other potentials. In reality, the induced electric field 
can be different if the boundary conditions that produce the rotating magnetic field are different. Though 
a rotating magnetic field is ubiquitous in circumstances ranging from an electric generator to a rotating 
neutron star jTD] , the treatment of the associated electric field can be a complicated problem that involves 
moving boundary conditions. Let us note however, that ^Bn{et) x r can be realized in a simple physical 
situation where the boundary is not moving at all. Consider two fixed and mutually perpendicular solenoids 
whose symmetry axes meet at the origin. The magnetic fields in the two solenoids have the same magnitude 
and change sinusoidally with 7r/2 phase difference. This produces a total magnetic field inside the common 
region that rotates about a third perpendicular axis through the origin. The induced electric field due to the 
change of magnetic field in each solenoid alone can be found by Faraday's law and by symmetry. The total 
induced electric field inside the intersection of the two solenoids can then be found to be equal to ^Bn{et) x r. 
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Adding a uniform and constant magnetic field along the third axis, we may obtain a uniform magnetic field 
that keeps at an arbitrary fixed angle with and rotates about the third axis whose vector potential is still 
given by ^Bn{et) x r. 

The induced electric field does not have to be ^Bh{€t) x r. Imagine for instance that the cross section 
of one of the solenoids is not a circle but instead a square, then symmetry argument no longer applies in 
calculating the electric field. One still produces the same uniform time dependent magnetic field, but the 
induced electric field is different and it is in general nonlinear with respect to the coordinates. So the vector 
potential is in general nonlinear. These different vector potentials should have different physical consequences 
but these will not be discussed further in this paper. 



3 The rotation transformation 

The time evolution operator U (t) satisfies 

iU{t)=HU{t), U{0)=I. (8) 

Here, we choose units in which h = c = 1 and we shall assume that repeated indices are summed over. 

Our starting point is to use the parallel transport method to transform to a rotating reference frame 
that follows the direction of the magnetic field. This rotating frame has the same origin as the initial 
coordinate system but with a rotating basis (ei, 62,63) which is fixed by 

ei{et) = (n X n) X e,{et), (9) 

where i — 1,2,3. The initial conditions are 

ei(0) = n'(0), e2(0) = n(0)xn'(0), 63(0) = n(0). (10) 

In the above, the prime in n'(0) means the derivative with respect to arc length s at t = while ej(et) denotes 
the time derivative. ( We assume that ds/dt > for all t.) It can be verified easily that 63 (et) = n(et) and 
that the three vectors form an orthonormal basis for all t. Furthermore the vectors ei(ei) and e2(ei) satisfy 
the condition ei(et) • 62 (et) — 0. Thus, when viewed as tangent vectors on the unit 2-sphere formed by 
all possible n's, ei(et) and 62 (et) have the geometric meaning of being parallel transported along the curve 
n{et) on the 2-sphere. Now accompanying such a transformation of the coordinate system we have 

U(t) ^ R{et)Ui{t) (11) 

where Ui describes the quantum time evolution in the rotating reference frame, and 

ft rniet) 
i?(et) = Texp ( - i / (n X h) • J dr) = Pexp (- i / (nxdn)-J). (12) 

Jo Jn{0) 

R{et) is a rotation operator well-known to be associated with the "solid angle Berry phase" [I]. It is 
characterized by the fundamental property that although it is only written as a time ordered exponential, it 
is in fact an explicit operator for a finite dimensional irreducible representation of the angular momentum 
and that if \m) is an eigenstate of n(0) • J with n(0) • 3\m) = m|m) and n{et) returns to n(0) at time T, 
then R{eT)\m) = exp{—imil.)\m) , where ft is the solid angle that the closed path traversed by n subtends 
at the origin of the parameter space. (For a review of this property, see section 2.1 of Ref [12].) 
The operator R has the property 

R-^{€t)e,{€t) -vRiet) ^ e, (0) • v ee , i = l,2,3, (13) 

where v is any vector operator. This property can be proven by first noting that it clearly holds for t = 0. 
Now it suffices to prove that the left hand side is constant, or its time derivative is zero. This is easily verified 
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by making use of the equations R = — i(n x h) • 3R, R~^ = R~^{t)i{n x n) • J, = (n x n) x e^, and the 
formula [a • J, b • v] = i(a x b) • v which holds for a vector operator v. From this property it is clear that 
Ui~^ViUiei{et) = J7"^fi[/ei(0). This is consistent with the fact that while the time evolution in the inertial 
reference frame is given by U{t), seen from the rotating frame, it is described by Ui. 
Using the above property, we can Taylor expand V{r ■ n{et) — L) and it follows that 

R-^{€f)V{Y-n{et)- L)R{et) = V{x-i-L). (14) 

Now we make the substitution U = RUi in the Schrodinger equation for U. It is straightforward to obtain 
the following equation for Ui 

iiliit) = HiUi{t), Ui{0) = I, (15) 
where Hi is the Hamiltonian in the rotating frame: 

= ^(p - eA(r))' + a2{et)Ji - ai(et) J2 + V{x3 - L) (16) 

with 

A(r) = A(r, 0) = Sn(0) x r/2, (17) 

and 

a^{et) = h{et) ■ e^(et), /x = 1, 2, J, = eijkXjPk- (18) 

The magnetic field is now in a fixed direction in the rotating frame. Other rotation operators can achieve 
the same goal but they introduce an additional term in Hi that is proportional to J3. The operator R leads 
to the simplest form of Hi. 



4 General factorization of U 

In this section our purpose is to find a factorization of Ui that has a simple form in the strong confinement 
limit. Suppose the confinement is strong, then motion of the particle in the direction of the magnetic field is 
dominated by the potential V{x3 — L). Also, the wave packet should be close to the x^ = L plane. However, 
directly extracting the operator exp[—i{pl/2m + V{x3 — L))t] from Ui does not lead to a useful factorization. 
This is because the term a2{et)Ji — ai{et)J2 in Hi contains both X3 and ps, and if we assign X3 = L in 
the limit of strong confinement, then is completely uncertain according to the Hcisenbcrg uncertainty 
principle. The issue here is that ps is not the kinematical momentum in the rotating frame in the current 
gauge. It is only after a gauge transformation that it is. This is the key observation that leads to the 
factorization of Ui. 

First let us rewrite Hi in the following form 

where Ki, K2 and K3 are the components of the kinematical momentum in the rotating frame, and Vc is 
the potential for the centripetal force: 

K^, = p^- eA^ - ma^{et)xs, fJ-=l,'2, (20) 
K3 = P3 + mai{et)xi + ma2{et)x2, (21) 
= ~h^{t)xl-^[ai{t)xi+a2{t)x2]\ (22) 

where A^, ji = 1,2 are given by A^ = A(r) • e^, or 

B B 

Ai = --^X2, A2 = —xi. (23) 
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In the rotating frame, in addition to the centripetal force, there are also the Coriolis and the Euler forces. 
These are produced by the terms —rna^ (et)a;3 and mai{et)xi + ma2{et)x2 in and K3. The term 
— eA^a^xa in Hi is responsible for the force due to the induced electric field caused by the rotation of 
the magnetic field and the Lorentz force arising from the rotation of the frame {ei{et),i = 1,2,3), i.e., the 
charged particle has the additional velocity (n x n) x r relative to the inertial frame {ei{0),i = 1,2,3) and 
this has a contribution to the Lorentz force. 

Our main step in elucidating the structure of Ui is to perform a gauge transformation as follows. Let |V'(0)) 
be an initial quantum state whose time evolution in the rotating frame is given by \il>i{t)) = Ui{t)\tp{0)) . 
Then i|^i(t)) = Hi\tpi{t)). In the rotating frame, we consider the gauge transformation 

\Mt))=9i^t)\Mt)), (24) 

where 

g{et) = exp[-im(ai(et)a;i + a2{et)x2){x3 - 2L)]. (25) 
In the gauge corresponding to |V'o(i))' the Schrodinger equation is 

i\Mt))=Ho\Mt)), (26) 

where 

Ho = g-\et)Hig{t) - ig-\et)g{et). (27) 
If the time evolution operator corresponding to Ho is Uo, then we have 

Ui=g{tt)Uog-\Q), iUo{t) = HoUo{t), Uo{0) = I. (28) 

Using the formula exp{—B)Aexp{B) = A+[A, B] with the condition that [A, B] commutes with A and B, 
we have 

^° " 2^"^ ^ 2^^^ ~ eA^a^xa + K - ig-\et)g{et) + -^pl + V{xs - L) (29) 

where 

= - eAfj, - 2ma^{et){x3 - L), 11= 1,2. (30) 

This implies that in the gauge corresponding to |i/'o(i))i the kinematical momentum is (IIi, 112,^3). 

The Hamiltonian Ho corresponding to this new gauge is unique in the sense that except i^p\, the rest 
of the terms in it do not depend on ps, and also Hi and 1X2 are in the simplest form in the limit of strong 
confinement. In fact, Hq can be rewritten as 

Ho = Hid + H2d + H^g , (31) 

where 

Hu = ^pI + V{X3 -L)- ^i^h^, (32) 

H2d = ^[{Pi - e^i)^ + {P2 - eA2f] - eLa^{et)A^ 

+ mL[ai{et)xi + a2{it)x2] - ^[ai{et)xi + a2{€t)x2]^, (33) 

H^o = [-^Oi^{et){pij, - eAij) - eA^a^(et) - md^(et)a;^]^o - mLi?^o + ^n^Co^, (34) 

with 

^0 = a;3 - L. (35) 
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If we now extract the operator 

Uu = exp I Hudt' = exp + V{x3 - L))t - [ ^L'^n^dt'] (36) 

Jo 2rn Jo 2 

from the operator Uq, such that 

Uo = Uuii, (37) 
iii = m, (38) 

S) = H2d + H^, (40) 



then it satisfies 

where = i?2d + U^J^H^gUid- Define 

then can be written as 
where 



= [-2a^(p^ - eA^) - ea^^A^ - md^(e<)a;^]^(f) - mLh^^{t) + ^n^^^t). (41) 
We see that ffj — >■ if ^ — 0. Now we can write il as 

il = U2dUi (42) 



iU2d = H2dU2d, U2d{0)=I, (43) 



where U2d is determined by 
and i7j is determined by 

iU^ = U^^'H^U2dUs, U^{0)=I. (44) 
Combining formulas, we obtain the factorization of the time evolution operator U{t) as follows: 

U{t) = R{et)Ui{t) = R{et)g{et)Uu{t)U2d{t)US)g-\Q). (45) 

This factorization holds generally, independent of the the specific form of the confinement potential and 
the rate at which n(ef) varies in time. However, if the confinement is strong enough, then ?7^^ff{[/2d ~ 
and ~ I. In this approximation, motion along the magnetic field direction and motion perpendicular to 
the magnetic field direction are described by two mutually commuting operators Uid and U2d- This is the 
main result of the paper. Also, based on the expression of in terms of ^{t), corrections to the strong 
confinement approximation may be obtained once U2d is known. 

Of relevance here is that X3 — L and C/j are operators, so the meaning of their closeness to the zero 
operator and the identity operator should be carefully analyzed. This is done in the following section where 
t/t is studied based on the expression of which has just been derived. Also, in the strong confinement 
limit — > /, U2d still corresponds to a general time evolution in which the rotation can be non-adiabatic. 
Only in the further limit 1/e m/eB, does it correspond to the adiabatic time evolution. The difference 
between the strong confinement condition and the adiabatic condition is an interesting point that is further 
examined in the following section. 

The main step in the factorization of Ui is the gauge transformation \ijji{t)) = g{et)\il)o{t)) , under which 
we switch from Ui to Uq through the relation Ui = g{et)Uog~^{0). In the gauge corresponding to Uq, the 
factorization is the easiest to derive: Uq = Uid{t)U2d{t)U^{t) . In the gauge corresponding to Ui, however, it 
is more natural to write the result as 

Ui{t) = Wu{t)W2d{t)W^it)g{et)g-\0), (46) 

where Wid = 9{(-i)Uid{t)g~^{et), and so on, still having the same physical contents and g{et)g~^{Q) serves 
to update the gauge of an initial state at t = to its gauge at time t. 
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5 and the strong confinement limit 



The operator ^(i) = U^^{t){xT, — L)Uid(t) (or, for the same purpose, ^0 = ^3 — L) is naturally associated 
with the strength of the confining potential. From the physical point of view, when the strength of the 
confinement potential increases, the energy gaps between different energy levels of the Hamiltonian Hid = 
2^p| + y(.X3 —L) also increase. For a given cigcnstatc |<I>„) of Hid, if {^n\x3 — L|<I>„) = 0, as is the case if the 
potential is symmetric with respect to the equilibrium position. thc;n the smallness of ^{t) can be measured 
by the quantity 

1^1 = Vm)) = v{^n\ixs-m^„). (47) 

If the minimum energy gap of Hid is A, then the strong confinement limit can be understood as cither |^| — > 
or A — >■ cxD. The two simplest examples are the harmonic oscillator potential and the infinite potential well. 
These correspond to linear and nonlinear Heisenberg equations of motion, respectively. For both cases, we 
have — c(n)/A, where c(n) depends only on the energy level n of Hid- Here we have three time scales, 
h/A, 1/e, and m/eB, which provide the ratios needed for a perturbative treatment of t/j. The condition 
that the confinement is strong means that ?l/A is much smaller than the other two time scales. 

Wc note that the factorization of U{t) here is not the same as ?7 = GDUc in that it reduces to a simple 
form in the strong confinement limit |^| — )■ rather than the limit e — > in which the rotation is infinitely 
slow. In the present problem, this factorization based on confinement strength is more general and more 
essential than a factorization based on rotation speed. This is because although H^ is proportional to both 
^{t) and e, the condition |^| — >• corresponds to [/{—>■ / but the condition e — >■ does not. To see this, let 
us examine the structure of U2dH^U2d, which determines U^. From the expression of H^, we see that the 
crucial term in U2^H^U2d is, 

U:^^[-2a^{et){p^ - eA^,) - ea^,{et)A^, - ma^{et)x^,]U2di{t). (48) 

This term causes the coupling between the ID and 2D motions. It is proportional both to ^(t) and e. It might 
seem that either the strong confinement limit A — > 00 or the limit e — > is sufficient in siippressing the effects 
caused by this term. But this is not so. The limit e — >■ is not enough. This is because the time evolution 
happens in the interval [0, 1/e], and terms of the order e may accumulate to give finite contributions. In 
the case that these e order terms are coupled to oscillating periodic functions with much shorter periods 
compared with 1/e, they are averaged out and do not accumulate. This is the mechanism behind the proof 
of the adiabatic theorem. Let us examine the term — 2Q!pJ72^^(et)(p^ — eA^)U2d£,{t) in the above expression. 
If e ^ 0, then U~J-{et){p^ - eA^,)U2d should recover the behavior of — eA^ as the Heisenberg operator in 
the usual Landau problem (see the next section), so it contains oscillating factors exp(±za;t). However, ^(t) 
is oscillating too and, in the simplest case of the harmonic oscillator, contains factors exp(±iAt//i). Then 
if Lj matches A//(,, such an averaging mechanism does not work. So e — > does not imply J/^ — > /. The 
strong confinement condition A — >• 00 is different because increasing A does not increase the duration of the 
time evolution. Thus, increasing A or, equivalently, decreasing |^| against fixed values of e and oj can make 
sure that the effects caused by U2dH^U2d are small for an initial state with fixed quantum numbers. So 
we deduce that the strong confinement condition is a more succinct statement here in considering the limit 

— )• I. Furthermore, in addition to shrinking |^|, increasing A against fixed values of e and ui automatically 
provides the fast periodic phase factors needed in averaging out the terms of the order e. This means that 

may in fact approach the identity operator faster than expected from the smallness of |^| alone. This 
may be useful in a more refined analysis on the difference between and the identity operator. 

Let us observe that if the potential is not close to the strong confinement limit, i.e., if the off-diagonal 
matrix elements of X3{t) — L are not small (i.e., the time scale h/A is not sufficiently smaller than 1/e and 
oj-i = m/(|e|B)), then a perturbation approach may become ineffective and the study of U^{t) may require 
different methods. It is clear that U^{t) does not in general commute with U^^ {t)x3Uid{t) because of its 
time dependence. Then U^(t) necessarily contributes to the dynamics both for motion along the magnetic 
field and motion perpendicular to it, with Uid and U2d not necessarily the major contributions. Also, U(^{t) 
is generated by the time-dependent Hamiltonian U2^H(U2d that is neither linear nor quadratic, because £,{t) 
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is in general a nonlinear function of X3 and p^. So the operator U^{t) contains the major complexity of the 
problem if the confinement is not strong. Nevertheless, is unique in the sense that it has concentrated the 
major complexities of the problem and that it approaches the identity operator in the strong confinement 
limit. For the general situation away from such a limit, perhaps this operator may contain useful structural 
information about the time evolution that can be extracted by methods beyond perturbation theory. One 
sees that U2d dictates a linear time evolution of the variables — eA^ and (/i = 1,2), so that the 
Hamiltonian U2^H^U2d represents a simple and nontrivial example of time-dependent coupling between 
linear and nonlinear motions. 



6 U2d and non-adiabatic effects 

Under the strong confinement limit, = I, the time evolution becomes 

U{t) = R{et)giet)Uidit)U2d{t)g-\0). (49) 

It is clear that in this limit, motion perpendicular to the magnetic field direction is determined by U2d{t). 
Although e does not have to be small compared with w, our focus here is to study U2d{i) by a perturbation 
method which is most effective for small e. 

Unlike -ff^, the Hamiltonian H2d is quadratic and corresponds to a linear system. For this kind of 
systems, there has been extensive work focusing on the construction of quadratic invariants and the quantum 
propagator (See, for instance, [T71 [SJ [TBI [13 ) ■ However, due to the existence of the term —^[ai{et)xi + 
a2(et)a;2]^ in H2d, the Heisenberg equations of motion are not solvable by quadrature for a general time- 
dependent function n(et) [20j . So in the present situation the quadratic invariants and the propagator cannot 
be constructed explicitly. 

Our approach here is to factorize from U2d a major contribution in terms of an operator that is explicitly 
constructed. Effects due to the existence of small terms including —^^[al{tt)xl + a2{et)x2Y are then studied 
by a perturbation method. In this way, U2d can be studied explicitly to any order in e. This approach also 
directly gives the Berry phase information contained in U2d- 



6.1 Factorization of U2d and the adiabatic limit 

The Hamiltonian H2d is of the form 

1 



H2d = ^[{Pi - eA^f + {p2 - eA2f] - eLa^{et)A^ + S{e^), (50) 

where S{e^) contains terms that are of the order e^: 

S{e^) = mLai^{et)Xfj, - y [ai(et)xi -|- a2iet)x2]'^- (51) 

Since the system is linear, effects caused by S'(e^) for the interval [0, 1/e] is at most of the order e, which 
goes to zero if e — )• 0. However, the term —eLa^{et)A^, which is of the order e, may accumulate over [0, 1/e] 
to have a finite effect. Denote 

H2d = ^[{Pi - eAif + {p2 - eA2f] - eLa^{et)A^, (52) 

U2d{t)^TcM-i f H2d{T)dT}. (53) 
JQ 

Then we have 

U2d = U2dU, (54) 
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where 



iJ, = Texp{-i/ U-^'Sie')U2ddT}. (55) 

JQ 

We first state a result about U2d{t)- Following Brown and Zak [131 [Ml [15] , we let 

TT^i = - eAf^, Vf^=Pfi + (iAf_,. (56) 

They satisfy 

[7rp,?7i.] = 0, [771,772] = -ieB, [7ri,7r2] = ieB. (57) 
The operators 77^, ^ — 1,2 generate magnetic translations: in addition to implementing 

exp(i7yidi + 17^2^2)2;^ exp(-i7;ic?i - 7772^2) = a;^ + d^, (58) 

just like ordinary translation operators, they commute with tt^. Assume e < 0, say it is the charge of the 
electron. Then the cyclotron frequency is a; = —eB/m. Using tt^ and 77^, one may define 

(7ri-77r2), 6 = (771 + 77/2). (59) 



Then a and b satisfy 

[a,at] = l, [b,b'']^l. (60) 
The usual Landau Hamiltonian can be written as 



— {ttI +Trl) = nuj{aa'< + 1/2). (61) 



In [16) , we studied the Landau problem with a time-dependent, spatially uniform electric field. The result 
in [12] directly applies here because —eLa^{et)A^^ is the same as a linear electric field potential. We have 

U2d = M{et)e-'"^'U,{t), (62) 

where M{et) is a path-ordered magnetic translation: 

M{et) = cxp[-i77idi(ei) - imd2iet)y'^'-'*\ (63) 

with dfj_{et) and I3{et) given by 

L /■* L 
df^i^t)^-- a^{eT)dT^-- • dn, (64) 

^ Jo ^ Jn(0) 

/3(rt) = = _eB( / di(er)( - ^:a2(eT))dT - idi(et)d2(ei))- (65) 

We see that Sd is the area enclosed by the path traversed by {di{et), di{et)) in the {di,d2) plane and the 
straight line pointing from {di{et), di{et)) to (0,0). Furthermore, it is clear that {di{et), di{et)) is fixed by 
the path of n{et) on the 2-sphere, so M{et) is determined by the path of n{et) on the 2-sphere. Therefore, 
M{et) is a geometric operator just as R{et) is. It corresponds to physically displacing a wave packet in the 
rotating frame through a path determined by {di{et),di{et)). The total length of the path is times the 
length of the path of n(ei) on the unit 2-sphere. The phase factor e*^ is quantum mechanical in nature and 
arises from the fact that the magnetic translation group is noncommutative. The operator J(t) is fixed by 
another path-ordered operator generated by tti and tt2, or equivalently, by a and a^: 

i7,(t) = e''(*)"-'^'W''^*^(*\ (66) 
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where 

5{t)^^j\ier)e-^-^dr, (67) 

j{t) = i f {S*SdT - SS*dT) = 4Ss. (68) 
Jo 

and 

h = l/V-2eB, a = ai+ia2. (69) 

Here Sg is the area enclosed by the path traversed by S{t) in the complex (5-plane and the straight line 
connecting the end and initial points of the path. But unlike d^, S is not determined by the path of n(et) 
alone. It also depends on how the path is traversed in time. Moreover, because of the oscillating factor e^'"^ 
in the integral, the integral goes to zero as S/u. So S should be considered as of the order e, and Ss is of the 
order e^. Therefore, Ue(t) — >■ / when e — !■ 0. In conclusion, U2d has a general expression 

U2d{t) = M{et)e~'"^'U,it)U,{t), (70) 

where Ue{t) -> /, Ue{t) — > / in the limit e — > 0. 

From this we see that the adiabatic limit, which assumes e — in addition to /, implies that 

U{t) = R{et)M{€t)e-'"^\ (71) 

Under the assumption of the harmonic oscillator potential V — iA:(r-n(ei) — L)^, this limit was derived before 
by analyzing the limiting behavior of the solution of the Heisenberg equations of motion in the rotating frame. 
However, the gauge transformation (7(e), which is a crucial ingredient in deriving the general factorization of 
U{t), was overlooked because g{e) I when e — ^ 0. 

In the adiabatic limit and for a cyclic variation of n(e<) with n(eri) = n(l) = n(0), the Berry phase 
factor is essentially i?(l)M(l). For different paths of n(ei), these i?(l)M(l)'s form a non-Abelian group, so 
the Berry phase factor is non-Abelian f2P . 



6.2 Non-adiabatic effects 

Among the consequences of the factorization of U{t), the operators and UJJe can be used in calculating 
corrections to the limiting cases. Here we study non-adiabatic effects under the strong confinement assump- 
tion that I. Then all non-adiabatic effects are in UfUe- Since is explicitly known, our focus here is 
the study of . 

Because S {e'^)U2d is of the order e^, [/^ — / is at most of the order e for i g [0, 1/e]. We have 

U,^I~i f U2^'S{e^)U2ddT + 0(e2). (72) 







Furthermore, because = I + 0(e), we have 

U, = I-i f e'"^^M-\eT)Sie^)M{eT)e-'"^^dT + Oie^). (73) 



Since in S'(e^) can be expressed in terms of a,a\b, and b\ the integral can be calculated by using the 
relation M{et)-'^Xf,M{et) ^ + d,, , Ug'^{t)hUB{t) = 6, and Ug'^{t)aUB{t) = ae"*"*. The result is 

[/, = (!- ici)I + C26 - c^6t + C366 - clb^b^ + €46^6 + cga^a + ©(e^), (74) 

where 

ci = [mLafj^d^ - —{aidi + a2d2 f - n^mlg]dT, (75) 
Jo ^ 
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/■* 1 

C2 = / -rnlBlLa* - {aidi + a2d2)a*]dT, (76) 
Jo ^ 

t 1 pt 

. 2 *2 



C3 — —i / —mlB a* dr, — C5 — i mlB n dr. (77) 
Jo 2 Jo 

In the above, the variables a^, a*, d^ and h under the integral sign are functions of t. In deriving this 

resuh, terms proportional to integrals such as a^(er)e~^*"'^rfr should be seen as of order because of the 

averaging effect of the oscillating factors. They are included in O(e^). The integrals in ci to C5 should be 

seen as of the order e because effects of the order in these expressions accumulate over the time interval 

[0, 1/e] and give effects of the order e. 

From the expression of U^{t), we see that the term S{e^) does not cause inter Landau level transitions 

up to the order e. Instead, the ci term represents a common phase correction for all wave functions, while 

csa^a represents a phase correction that depends on the Landau energy level. The other terms represent 

translations (c2 and C2) and distortions (03,03 and C4) of a wave packet without changing its energy, which 

are corrections to the non-Abelian part of the geometric phase. The inter-Landau level transitions caused by 

the term —eLa^{et)Af^ as represented by C/e is equivalent to the influence of a uniform electric field which 

has been analyzed in detail in |16) . 



7 Conclusions 

In this paper we have studied in detail the Landau electron problem in a rotating magnetic field and a 
rotating potential. The Landau energy levels of the instantaneous Hamiltonian are infinitely degenerate and 
the potential is of a general form. The key result of the paper is a natural factorization of the time evolution 
operator by making use of the rotation operator technique and by a careful analysis of the structure of the 
Hamiltonian in the rotating frame. The factorization is from a different perspective than the one suggested by 
the quantum adiabatic theorem and includes the latter as a special case. The detailed structure regarding 
the coupling between motion in the magnetic field direction and motion perpendicular to it is carefully 
analyzed. It is pointed out that it represents a simple and nontrivial example of the coupling between linear 
and nonlinear motions. 

In the strong confinement limit, the problem is reduced to a two-dimensional one. We presented a 
different method so that effects due to the quadratic potential can be analyzed perturbatively. The time 
evolution operator for this two-dimensional motion is completely determined to include all e order terms. 

Since the Landau electron problem and the rotating magnetic field are often encountered in different 
areas in physics, the method and results presented here may find their use in the study of various physical 
situations associated with applications. 



References 

[1] M. V.Berry, Proc. R. Soc. A392, 45-57 (1984). 

[2] Geometric Phases in Physics, edited by A Shapere and F Wilczek, World Scientific (1989). 

[3] For an introduction to some of the applications of the geometric phase, see A. Bohm, A. Mostafazadeh, 
H. Koizumi, Q. Niu, and J. Zwanziger, Tiie Geometric Phase in Quantum Systems, Springer- Verlag, 
Berlin (2003). 

[4] A. Messiah, Quantum Mechanics (North Holland, Amsterdam, 1970), Vol. 2. 

[5] J. E. Avron, R. Seller, L. G. Yaffe, Commun. Math. Phys. 110 33 (1987) (Erratum: Commun. Math. 
Phys. 153, 649-650 (1993)). 

[6] V.V. Dodonov, LA. Malkin and V.I. Manko, Physica 59 241 (1972). 



12 



J. Chee, Phys. Lett. A 275 473 (2000). 

G. Date and P. P. Divakaran, Ann. Phys. (N.Y.), 309, 421 (2004). 
M. H. Al-Hashimi, U.-J. Wiese, Ann. Phys. 324 343 (2009). 

M. Bocquet, S. Bonazzola, E. Gourgoulhon, J. Novak, Astron. Astrophys. 301, 757 (1995) 
M.V. Berry, Proc. R. Soc. A414, 31-46 (1987). 
J. Chee, J. Phys. A: Math. Theor. 43 345301 (2010). 
E. Brown, Phys. Rev. 133, A1038 (1964). 
J. Zak, Phys. Rev. 134, A1602 (1964). 
J. Zak, Phys. Rev. 134, A1607 (1964). 
J. Chee, Ann. Phys. 324 97 (2009). 

H. R. Lewis and W. B. Riesenfeld, J. Math. Phys. 10 1458 (1969). 

M. Combescure and D. Robert, Quadratic Quantum Hamihonians Revisited, arXiv:math-ph/0509027^ 1, 
(2005). 

G. Fiore and L. Gouba, J. Math. Phys. 52, 103509 (2011). 

R.L. Bryant, An Introduction to Lie Groups and Symplectic Geometry, Lectures delivered at the Re- 
gional Geometry Institute in Park City, Utah, 24 June-20 July 1991, pages 51-52 and problem 10 on 
page 59. 

[21] F.Wilczek and A. Zee, Phys. Rev. Lett. 52 2111 (1984). 



13 



